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Preface 


Literature  on  the  many  methods  of  manual  orbit  deter¬ 
mination  Is  extensive  as  one  might  expect.  However,  in 
the  course  of  my  research  for  simple  computational  methods 
and  devices,  I  noted  that  the  graphical  device  known  as  a 
velocity  hodograph  had  virtually  disappeared  from  recent 
literature.  I  felt  and  still  feel  that  the  velocity  hodo¬ 
graph  offers  interesting  possibilities  for  solving  orbit 
determination  problems.  To  a  large  extent,  this  thesis  Is 
devoted  to  practical  applications  of  the  velocity  hodo¬ 
graph.  Since  this  graphical  device  incorporates  the  basic 
concepts  of  orbital  mechanics,  I  claim  no  truly  original 
Ideas . 

I  would  like  to  take  this  opportunity  to  express  my 
appreciation  to  Professor  Peter  Blelkowlcz  for  leading  me 
through  several  courses  in  astrodynamlos,  suggesting  this 
topic,  and  finally,  but  certainly  not  least,  for  his 
guidance  and  criticism  during  the  writing  of  this  thesis. 

I  also  wish  to  acknowledge  my  gratitude  to  Captain  Tony 
Barth  of  the  USAF  Avionics  Laboratory  for  providing  me 
with  many  research  documents.  My  special  thanks  to  my 
family  for  their  understanding  and  encouragement  during 
the  preparation  of  this  report,  with  particular  expres¬ 
sions  of  appreciation  to  my  wife,  Evelyne,  for  her 
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Abstract 

Simple  non-electrical  methods  and  devices  for  orbit 
determination  on  the  basis  of  explicit  navigation  are 
examined.  Principal  measuring  devices  discussed  are  the 
sextant,  stadimetrio  instruments,  and  almucanters.  Com¬ 
putational  methods  for  solving  the  geometric  orbital 
parameters  are  focused  on  adaptations  of  the  velocity  hodo- 
graph.  Discussion  Includes  sample  problems  and  a  general 
analysis  of  inherent  errors. 
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METHODS  AND  DESIGN 
OF  SIMPLE  DEVICES  FOR  ON-BOARD 
AND  GROUND-BOUND  ORBIT  DETERMINATION 


I.  Introduction 


Background 

After  a  decade  of  successful  space  exploration  In  the 
vicinity  of  the  earth  which  was  climaxed  by  the  Apollo 
landing  on  the  Moon,  man  continues  to  strive  to  penetrate 
deeper  into  space.  Of  particular  Interest  among  possible 
future  manned  missions  are  fly-bys  and  landings  on  other 
near  and  distant  planets  of  the  solar  system.  Due  to  the 
great  distances  Involved,  Earth-based  radio  navigation 
capabilities  will  be  severely  limited  if  for  no  other  rea¬ 
son  than  time-delayed  communication.  Intuitively,  this 
situation  leads  to  an  emphasis  on  developing  on-board 
navigation  procedures  and  measuring  devices  which  are 
totally  Independent  and  self-sufficient  with  respect  to 
collection  of  data  (the  observables)  and  orbit  computa¬ 
tions.  In  addition,  missions  landing  on  some  planets  may 
be  faced  with  determining  the  orbits  of  probes  or  command 
modules  from  observations  made  on  the  ground.  The  method 
which  presents  the  greatest  reliability,  although  not 
necessarily  the  greatest  accuraoy,  is  most  often  the  one 
divorced  from  sophistication.  Since  electronic  equipment 
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failure  is  a  distinct  possibility,  simple  non-electrical 
devices  augmented  by  manual  calculation  schemes  which  are 
most  general  in  applicability  should  be  considered. 

Problem 

This  report  is  intended  to  present  a  discussion  of 
methods  and  design  of  simple  devices  for  on-board  naviga¬ 
tion  and  ground  level  orbit  determination,  even  In  the 
event  of  electronic  equipment  failure.  "On-board"  naviga¬ 
tion  is  defined  as  space  navigation,  the  art  of  describing 
the  spacecraft  state  from  observables  as  seer,  from  within 
the  spacecraft,  without  resort  to  outside  communication. 
"Ground  based  orbit  determination"  is  defined  as  the  art 
of  describirg  a  spacecraft  state  from  observables  as  seen 
from  the  ground.  The  term  "state"  is  defined  as  the  posi¬ 
tion  and  velocity  relative  to  some  reference  coordinate 
frame  for  a  particular  point  in  time.  Classically,  the 
state  is  expressed  by  a  set  of  numerical  constants  called 
orbital  parameters  when  the  motion  of  the  vehicle  is  in  a 
free-fall  trajectory. 

Scope 

Prev‘ jUs  theses  dealing  with  manual  navigation  were 
written  by  Captains  Schehr  and  Smith  ("Manual  Astronaut 
Navigation!  Apollo  Mission  Application")  and  by  Major 
Horrigan  and  Captain  Walsh  ("Manual  Astronaut  Navigation") 
under  the  sponsorship  of  Professor  Bielkowicz.  These 
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gentlemen  developed  the  Delta-H  method  which  was  Initially 
applicable  to  small  range,  low-eccentrlclty  orbits  and  then 
extended  to  highly  eccentric  elliptical  orbits.  In  addi¬ 
tion,  these  theses,  combined  together,  form  a  wonderful 
synopsis  of  the  measurements  which  can  be  made  with  a 
hand-held  sextant  or  stadimeter  and  the  computations  which 
follow  to  specify  the  orbital  parameters. 

Rather  than  repeating  their  voluminous  Information  In 
new  or  disguised  verbage,  the  principle  objective  is  to 
Investigate  simple  non-electrlcal  devices  which  can  be 
used  to  collect  observable  data  and  the  promising  manual 
methods  of  calculation  using  the  velocity  hodograph. 

Assumptions 

The  primary  assumption  made  is  that  the  spacecraft 
follows  unperturbed,  pure  Keplerian  flight.  This  Implies 
that  i 

1.  The  trajectory  is  restricted  to  two-body 
motion  and  describes  conic  sections. 

2.  The  gravity  field  of  the  attracting  body  is 
uniform.  Inverse-square,  and  central. 

3.  The  mass  of  the  spacecraft  is  negligible 
compared  to  the  attracting  body  mass. 

4.  The  vehicle  is  always  consiiered  to  be 
located  In  one  clearly  defined  sp  e  of 
influence  such  that  travel  from  one  body  to 
another  Involves  patch  conic  solutions. 
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Other  assumptions  made  are  that  there  are  no  pre-cal- 
oulated  reference  trajectories  which  therefore  Implies 
relying  on  explicit  navigation  schemes,  and  that  the 
reader  is  acquainted  with  the  basics  of  orbital  mechanics. 

Organization 

Chapter  II  deals  with  the  design  of  simple  devices 
used  to  measure  the  observables  from  a  space  vehicle. 

Chapter  III  presents  a  theoretical  discussion  of  the 
velocity  hodograph. 

Chapter  IV  provides  examples  of  typical  problems  and 
their  solution. 

Chapter  V  deals  with  grov  id  level  orbit  determination 
and  associated  hardware. 
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o 

II.  Simple  Measuring  Devices 


Introduction 

Space  navigation  as  well  as  ground  level  orbit  deter¬ 
mination  are  both  an  Instrumentation  and  computational 
problem.  All  Instrumentation  and  computational  modes 
should  employ  the  same  navigation  concept.  The  concept 
selected  In  this  report  Is  the  "explicit"  navigation 
scheme.  Explicit  navigation  is  more  desirable  In  the 
sense  that  It  can  handle  any  situation  which  calls  for  a 
deviation  from  the  flight  plan  (Ref  12t428).  Although 
the  explicit  navigation  method  has  the  disadvantage  of  a 
large  computational  load.  It  has  the  primary  advantage  of 
providing  the  navigator  with  direct  data  on  spacecraft 
position  and  indirectly  the  velocity.  The  raw  data  that 
must  be  obtained  for  use  In  the  computations  are  position 
fixes.  The  position  fix  data  Includes  a  radial  line  of 
position  (LOP)  from  the  attracting  body  and  a  radial  dis¬ 
tance  on  the  LOP.  Three  position  fixes  permit  the  naviga¬ 
tor  to  fully  define  the  orbit.  Partial  fixes  which 
include  only  the  LOP  or  the  radial  distance  can  also 
establish  the  orbital  parameters  completely  or  partially 
depending  on  the  conditions  and  circumstances)  for  example, 
having  only  three  ranges  would  only  allow  computing  three 
orbital  parameters  (See  Chapter  III  for  more  details). 
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Figure  1.  The  Sextant 
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The  iype3  of  non-electrical  instruments  that  can  provide 
LOP  and  range  (radial  distance)  Information  are  all  opti¬ 
cal  instruments.  Fortunately,  the  sighting  and  aiming 
Involved  is  one  area  where  the  superior  capabilities  of 
the  human  eye  for  resolution  and  recognition  can  be  used 
advantageously.  Discussion  of  the  individual  Instruments 
now  follows. 

The  Sextant 

The  sextant  is  an  Instrument  fitted  with  moving 
mirrors  and  prisms.  The  basic  measurements  taken  are 
angles  between  the  line-of -sight  of  one  reference  body  to 
another,  or  to  measure  the  angle  subtended  by  the  diameter 
of  an  apparent  disc  such  as  the  sun,  moon,  or  planets. 

Oddly  enough,  the  ancient  mariner’s  sextant  (Fig.  1) 
is  readily  adaptable,  in  principle,  to  space  navigation 
and  orbit  determination  problems.  This  adaptability  was 
dramatically  demonstrated  <.uring  the  manned  space  flights 
of  Gemini  IV,  VII,  and  XII. 

The  required  accuracy  of  sightings  for  lunar  or 
planetary  missions  is  several  seconds  of  arc  (Ref  4t208). 
Current  aiarine  sextants  such  as  the  Havy  Mark  II  Mod  0 
weights  approximately  2,7  pounds.  The  vernier  readout 
least  count  is  0.1  minute  of  arc  (Ref  9*3) •  Although  it 
is  a  well  known  fact  that  simply  increasing  the  radius  of 
the  sextant  arc  will  correspondingly  increase  the  instru¬ 
ment  accuracy,  the  limitations  placed  on  weight,  size, 
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Figure  2.  Schematic  Diagram  of 
Optical  Hand-held 
Sextant 

(From  Ref  20  «4) 

restricted  field  of  view,  and  cramped  quarters  in  the 
current  space  program  makes  Instrument  accuracy  purely  an 
engineering  problem.  Thus,  the  simple  sextant  becomes  a 
sophisticated  piece  of  hardware  containing  such  items  as 
precision  stainless  steel  worm  and  worm  wheels  (35  mil¬ 
lionth  of  an  inch)  and  beryllium  structural  elements  (Ref 
14i198).  The  complexity  may  be  easily  visualized  by 
glancing  at  Figure  2. 

Studies  made  at  the  NASA  Ames  Research  Center, 


Moffett  Field,  California,  slu  ed  that  the  ability  of 
sextant  operators  to  make  accurate  navigation  measurements 
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was  critically  dependent  of  training.  Test  results  at  the 
research  center  also  Indicated  that  Increasing  the  sextant 
telescope  magnification  frcr:  2  to  20  power  decreased  the 
average  standard  deviation  from  22  seconds  of  arc  to  about 
5  arcsec.  The  results  of  Gemini  XII  showed  that  hand-held 
sextant  performance  was  excellent,  providing  a  standard 
deviation  measurement  error  of  less  than  +10  arcsec  (Ref 
19*655). 

A  photographic  sextant  that  can  make  conventional 
sextant  measurements  and  simultaneously  provide  positive 
means  for  target  star  identification  was  also  investigated 
at  the  Ames  Research  Center.  The  summary  of  the  study 
concluded  that  the  film  was  the  most  critical  factor  and 
that  currently  available  emulsions  lack  either  resolution, 
speed,  or  ease  of  development  (Ref  24«1), 

The  Stadlmeter 

The  stadlmeter  is  a  device  used  to  obtain  the  dis¬ 
tance  separating  the  observer  from  another  body  by  taking 
angular  measurements  of  that  body  and  provided  the  body’s 
diameter  is  known.  There  are  two  classes  of  systems i 
systems  viewing  a  limited  portion  of  the  target  and  sys¬ 
tems  viewing  the  entire  target. 

Development  of  the  stadlmeter  by  the  U.S.  Air  Force 
and  the  Kollsman  Instrument  Corporation  has  emphasized 
the  class  of  instruments  measuring  the  amount  of  curvauure 
evident  in  the  limited  field  of  view  available.  The 
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observed  curvature  can  be  translated  into  range  by  means 
of  precomputed  calibration  curves.  One  way  to  obtain  the 
curvature  is  to  measure  the  deviation  angle  between  two 
chords  connecting  three  equidistant  horizon  points  (A,  B, 
and  C).  A  pair  of  composite  prisms  creates  three  super¬ 
imposed  fields  of  view.  The  deviation  angle  is  measured 
when  the  three  arcs  of  the  horizon  intersect  (Fig.  3)» 
Another  method  consists  of  a  triple  split  image  field  with 
a  center  scanning  prism  (Fig.  4)  which  when  aligned  with 
the  outer  segments  1  and  3*  the  angle  0  can  be  read  di¬ 
rectly  from  the  stadlmeter. 

The  equation  providing  the  relationship  between  the 
measured  angle  0  and  the  altitude  h  as  geometrically  pre¬ 
sented  in  Figure  5  is 


1}  (2-1 )* 


where 

R  =  radius  from  body  center  to  horizon 
A  *  fixed  angle  between  the  side  lines  of  sight 
0  =  measured  angle  0  plus  a  correction  angle  J0 
The  pre -computed  calibration  curves  plot  h  versus  0»  how- 
ever,  from  Eq  (2-1)  it  is  evident  that  (h/R)  versus  0  can 

♦This  equation  is  derived  in  Ref  I8tll7» 
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INSTRUMENT  FIELD  OF  VIEW 


1.  2.  3. 


1 .  Stadimeter  pointed  at  horizon  so  that  the  images 
of  the  horizon  at  A  and  B  cross  horizon  at  C  out 
of  the  field  of  view. 

2.  Horizon  at  C  brought  into  field  of  view, 

3t  Horizon  at  C  aligned  upon  the  intersection  of  the 
horizons  at  A  and  B.  The  reticle  aids  in  this 
operation. 


Figure  3»  Horizon  Geometry  and  Alignment  of 
Images  in  Instrument  Field  of  View 
(From  Ref  25  * 8—9) 
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Figure  5»  Geometry  of  Stadimetric  Measurement 

(From  Ref  13*103) 


also  be  plotted.  Since  0  and  £9  are  both  dimensionless 
functions  of  R  and  h,  a  plot  of  (h/R)  versus  ©would  be 
more  general  since  it  would  apply  to  any  heavenly  body 
(See  Appendix  A  for  specific  relatioships) . 

The  stadimeters  Just  discussed  are  intended  for  use 
at  low  altitudes-  For  example,  the  USAF  Space  Stadimeter 
is  not  Intended  for  use  above  2500  NM  altitude.  For 
spacecraft  altitudes  above  2500  NM,  the  angular  diameter 
of  the  body  can  be  measured  (Fig.  6).  The  following 


Figure  6.  Stadlmetric  Angle  With  Respect 

to  a  Spherical  Body 


equation  is  used  to  determine  range 


d/a 

s  in  % 


R 

Sot  P/z 


(2-2) 


where 

d  =  diameter  of  body 

r  ngle  subtended  over  diameter 
The  second  class  of  stadimeters  views  the  entire 
target.  A  typical  method  is  one  proposes  by  L.M.  Vorob'ev 
(Ref  23«126)  which  consists  of  a  movable  circular  screen 
(Fig.  7)  whose  plane  is  perpendicular  to  the  direction  of 
the  celestial  body.  The  screen  of  diameter  2r  moves  per¬ 
pendicularly  to  the  axis  00,  to  a  ">sition  where  the  edges 
of  the  celestial  body  are  barely  visible.  The  subtended 
half  angle  from  Figure  7  is 
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Figure  ?.  Measurement  of  the  Angular 
Diameter  J&  of  a  Planet 
By  Means  of  a  Movable 
Circular  Screen  of  Radius  r 
(From  Ref  23*126) 


(2-3) 


Substituting  the  value  for  the  planet  radius  and  the 
angle  3) 2  Into  Eq  (2-2)  produces  the  range  value. 

Zenith  Viewing  Telesoope 

A  simple  optical  Instrument  called  a  "bazooka*'  or 
zenith  viewing  telescope  as  illustrated  in  Figure  8  can  bs 
used  to  establish  the  local  vertioal  of  the  spaceship 
above  the  surface  of  the  celestial  body  (Ref  l6i289)» 
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Figure  8.  Zenith  Viewing  Telescope 
(From  Ref  26*790) 

Theoretically,  this  method  amounts  to  star-matching.  An 
Identified  star  at  the  zenith  would  immediately  establish 
the  line  of  position  in  terms  of  right  ascension  and 
declination. 

The  Almuoantar 

The  almucantar  is  a  general  application  of  the 
"bazooka'*  telescope.  This  system  calls  for  establishing 
a  local  vertical  vector  which  is  the  symmetric  axis  for 
cones  whose  semi-angle  is  measured  from  the  zenith.  Thus 
a  simple  peep  sight  centering  element  placed  at  the  cone 
apex  and  concentric  circles  scribed  on  a  window  with  the 
zenith  at  the  center  would  theoretically  suffice  (Ref  10 1 
249) .  The  concentric  circles  represent  degrees  of  co-alti- 
tude.  Knowing  the  cc-altitude  of  two  stars  at  any  Instant 
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of  time  can  establish  a  LOP,  (assuming  that  the  navigator 
can  eye  ball  the  zenith  to  resolve  any  ambiguity). 

Angular  Accuracy 

In  reality,  arv  parameter  Is  measured  with  some  error. 
Besides  the  usual  Instrumentation  and  personnel  error,  such 
factors  as  the  aberration  of  light,  confusion  with  binary 
stars,  and  proper  motion  of  stars  can  contribute  to  erro¬ 
neous  measurements. 

Aberrational  displacement  due  to  the  spacecraft's 
velocity  Is  small  for  velocities  of  10  to  20  km/sec  but 
for  velocities  of  about  100  km/sec,  the  angular  displace¬ 
ment  may  approach  1  minute. 

The  distance  between  components  of  binary  stars  are 
relatively  large.  For  example,  the  components  of  the 
binary  star  In  Ursa  Major  a.e  situated  about  12'  apart 
(Ref  23il25).  The  best  solution  for  alleviating  this 
error  is  to  avoid  using  binary  stars. 

Proper  motion  of  stars  can  be  dismissed  on  the  assump¬ 
tion  that  the  star  catalog  In  use  is  current  and  that  the 
annual  proper  motion  of  stars  is  measured  In  mlcroradlans . 

The  effect  of  angular  measurement  error  increases  with 
the  distance  from  the  body.  For  example,  a  space  vehicle 
situated  approximately  80  million  miles  from  the  sun  and 
20  million  miles  from  each  of  two  planets  being  sighted, 
a  total  error  of  30  arosec  would  result  In  a  12,000  mile 
error  in  position,  In  a  direction  perpendicular  to  a  line 
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to  the  sun  and  alxmt  3,000  miles  In  directions  perpendic¬ 
ular  to  lines  connecting  the  vehicle  with  the  planets 
(Ref  809). 

Range  Accuracy 

From  Eq  (2-2)  it  is  evident  tnat  for  a  given  angular 
error  &  and  a  given  body  radius  R  that  the  range  error  is 
magnified  as  the  range  Increases. 

The  more  difficult  error  to  handle  is  the  one  con¬ 
cerning  the  body  radius.  Errors  in  body  radius  measure¬ 
ment  are  attributable  to  the  influence  of  the  planets  * 
atmospheres  and  their  nonsphericity.  To  a  great  extent 
the  former  reason  is  due  to  a  question  of  definition. 

For  example,  the  radius  of  the  Earth  for  stadlmetrlc  meas¬ 
urements  is  taken  as  the  alrglow  radius  which  is  14.9  NM 
above  the  average  earth  surface  radius.  Since,  as  previ¬ 
ously  mentioned,  this  error  Is  due  to  definition  It  will 
not  be  discussed  further. 

Looking  at  the  non-sphericity  problem,  the  assumption 
may  be  made  that  if  the  body  is  not  a  sphere  then  It  is  an 
ellipsoid  of  revolution  about  its  spin  axis.  Just  as  It 
Is  necessary  to  know  the  radius  of  a  spherical  body,  here 
the  oblateness  must  also  be  known j  therefore,  let  It  be 
assumed  that  the  following  physical  characteristics  of  the 
planet  are  common  knowledge 

a  =  semi-major  axis 
b  =  semi-minor  axis 
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e  =  planet  eccentricity 

Since  the  planet  is  like  a  spinning  top  or  gyro,  the 
intercept  of  the  celestial  sphere  by  the  axis  if;  fixed  in 
space  and  can  be  identified  or  at  least  nearly  approximated 
by  some  star.  For  *130510  In  Figure  9»  the  intercept  is 
labeled  as  Star  Q.  Angle  t  is  the  smallest  or  largest 
angle  measured  from  Star  Q  to  the  body  surface  In  the  XY 
plane  or  any  other  plane  containing  the  spin  axis.  Look¬ 
ing  at  the  ellipsoid  from  a  great  distance  the  body  ap¬ 
pears  as  a  circle  of  radius  a  when  f  «  180°  and  it  ap¬ 
pears  as  an  ellipse  with  a  semi-minor  axis  equal  to  b 
when  f -  90°.  In  fact,  from  any  intermediate  angle  Jf, 

the  body  seen  from  a  great  distance  such  as  at  point  A 
appears  as  an  ellipse  whose  semi-minor  axis  R  is 

ft  =  a  (l  -  fSsLnt'}'^  (2-4) 

Note  that  R  is  also  the  apparent  radius  of  the  body  meas¬ 
ured  in  the  XI  plane.  The  point  of  tangency  is  A*,  but  a 
perpendicular  constructed  from  the  line-of-slght  APA'  to 
the  body  center  is  R,  (the  apparent  radius  as  viewed  from 
point  A) . 

At  point  P  the  angles  If,  and  or  X,  and  <f>,  may  be 
measured  with  a  sextant  and  the  sliding  screen  stadlmeter. 
R,  and  Rt  can  be  calculated  from  Eq  (2-4).  Since  <f>  is  the 
difference  between  if,  and  suppose  that  if,  and  <f>  are 
the  angles  measured.  From  the  problem  geometry 
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where 


Y' 


R^k  = 


ft. 

Sin  «<, 


sln(4>  -<) 


(2-5) 


r  =  radial  distance  to  vehicle 

Bfc  =  body  surface  radius  on  the  line  joining 
the  body  center  and  spacecraft 

h  =  altitude 

°f,  =  subtended  radius  arc  for  the  apparent 
radius  B, 


Application  of  a  trigonometric  identity  and  rearrang¬ 
ing  Eq  (2-5)  yields 


( 


sin  <t>  cos <*,  -  Cos <t>  s in. °<l 
sin.oi, 


(2-6) 


therefore 


o< 

i 


(2-7) 


Having  solved  o<(,  let  ^/2  =  c<f  and  substitute  into  Eq 
(2-2)  to  obtain  the  range  r. 

The  angle  A.  which  corresponds  to  the  body  latitude 
may  be  found  from  the  following  equation 
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71  *  /SO*  -  (<+«<,)  (2-8) 

Knowing  A#  the  body  radius  Bj,  on  the  line  connecting  the 
body  center  with  the  spacecraft  is 

O  _  ^  ^ _  (2-9) 

**  smlA.  +  loz  cosZ a)/z 

and  the  altitude  can  be  immediately  obtained  from  Eq.  (2-5) • 
As  just  shown,  the  problem  of  non-sphericity  can  be 
handled  through  additional  calculations  so  that  oblateness 
need  not  introduce  range  errors.  Thus,  the  accuracy  in 
range  is  really  dependent  on  how  well  the  physical  dimen¬ 
sions  of  a  celestial  body  are  known. 

u 
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I.'.I.  The  Velocity  Hodograph 


Introduction 

Graphical  solutions  can  he  developed  advantageously 
In  the  solution  of  problems  In  space  flight.  Graphs  are 
pictures,  whieh,  as  Is  so  often  said,  are  worth  a  thousand 
words.  It  geometrically  portrays  numerical  results  that 
are  mathematically  difficult  to  grasp.  The  velocity  hodo- 
graph  defined  as  the  locus  of  velocity  vectors  Is  a  geo¬ 
metrical  picture  which  reduces  all  conic  sections  to  map¬ 
ped  circles. 

Mathematical  Derivation 

If  motion  of  the  spacecraft  Is  described  In  polar 
coordinates,  the  velocity  vector  can  be  defined  by  two 
orthogonal  components i  the  radial  velocity,  r,  and  the 
transverse  velocity,  r0  .  The  equations  used  to  develop 
the  hodograph  equation  are 

V2  =  rS  (Ve)2  (4  -  i-J  (3-1) 

kz  - y*.  [a.  (l  -  e*)j  (3-2) 

k  =  T’O  -  ConsW  (3-3) 
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where 


a  =  semi-major  axis 
e  *  eccentricity 

h  =  angular  momentum  per  unit  mass 

V  =  total  velocity 

r  =  radial  distance 
* 

Q  -  angular  velocity 
yM-  gravitational  constant 
Substituting  Eq  (3-3)  into  (3-2)  and  rearranging 


yields 


0.  "  k& 


(3-4) 


Eq  (3-4)  can  now  be  directly  substituted  into  Eq  (3-1)  to 


give 


f%  (r8)1  >>(-.? 


(3-5) 


which  when  rearranged  yields  Eq  (3-6)  below.  It  is  imme¬ 
diately  recognized  as  the  equation  of  a  circle  with  coor- 

« 

dinate  axes  r  and  r0  ,  It  should  be  noted  that  in  the 
proceeding  development  the  initial  equations,  Eq  (3-1)  and 
Eq  (3-2),  are  equations  for  circular  and  elliptical  sec¬ 
tions;  however,  the  same  resulting  equation  is  produced 


assuming  parabolic  or  hyperbolic  conic  sections. 


(3-6) 
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Eq  (3-6)  as  notated  represents  a  graph  with  a  "float 
ing"  center *.  It  is  desirable  to  fix  the  center  snch  that 
a  pre-plotted  graph  can  be  constructed  and  used  under  all 
circumstances .  This  can  be  accomplished  by  reducing  the 
equation  to  a  non-dimensional  form.  Thus,  dividing  Eq 
(3-6)  by  (>“/h f  yields 


(3-7) 


Letting 


y**- 


(3-8) 


and 


(3-9) 


it  is  evident  that  in  the  X,Y  coordinates  the  origin  for 
circles  of  radius  e  is  centered  at  the  point  (1,0)  as 
shown  in  figure  10. 


Discussion 

Many  interesting  relationships  can  be  obtained  by 
just  looking  at  X  and  Y,  for  example t 

,  .  .  Crdf.  _kZ 


x  =  — - 


Vs- 

•lc 


y4*"  v” 


=  /4ec  osQ  (3-10) 
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Figure  10.  Dimensionless  Polar 

Hodograph  of  Velocity 


and 


Y  =  ^  - 

i  ~ 


'e  .  ^X; 


=  e  sin.9 


(3-U) 


where  v„.  the  local  circular  velocity,  is Obviously 
many  more  rearrangements  are  possible.  From  Eq  (3-10)  it 
becomes  immediately  apparent  that  9  is  the  true  anomaly 
and  since  X  is  equal  to  (1+e)  when  9  is  zero  and  (r@)  is 
maximum,  it  follows  that  the  point  (l+e,0)  is  the  peri¬ 
focal  point.  By  the  same  type  of  argument  (l-e,0)  is  the 
apofocal  point;  however,  it  should  be  noted  from  Eq  (3-8) 
that  points  where  X  is  equal  to  or  less  than  zero  are  not 
physically  defined.  Furthermore,  it  should  be  noted  that 
the  X  axis  represents  the  apsidal  axis.  As  the  true  anom¬ 


aly  increases  from  perifocal  point,  r  is  always  initially 
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positive  so  from  Eq  (3-lD  i*  Is  deduced  that  6  is  meas¬ 
ured  in  the  counterclockwise  direction  since  Y  is  initial¬ 
ly  positive. 

Additionally,  it  is  evident  that  the  arctangent  of 
(Y/X)  is  the  same  as  (r/r0 )  so  that  a  line  Joining  point  0 
to  point  (X,Y)  forms  an  angle  with  the  X  axis  equal  to  the 

flight  path  angle. 

Time  from  peripoint  passage  can  be  indirectly  obtain¬ 
ed  by  adding  new  curves  to  the  graph.  Consider  the  time 
equations  (Ref  li4-5)  for  an  ellipse, 


I  + 


sin  9 
ecoiQ  1 


(3-12) 


for  a  hyperbola, 


At-* 


sin  9  _  +  i/e-l  tan  t 


i  +  ecos0  '  ^--4 Sh  tan- £ 


f 


V  (3-13) 


and  for  a  parabola 


(3-14) 
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Noting  that 


UM*,X 

\  ^  s*t 


Eq  (3-12),  (3-13) »  and  (3-1^)  can  be  written  as 

t(X  -  KrX^ 


(3-15) 


(3-16) 


where  K  =  f  ( e ,  0  ) ,  thus 


K  =  0 


for  0  =  0* 


•m 

K  =  Ttjl-e1)  2  for  0  =  180° 

By  preplotting  constant  lines  of  K  (Ref  3:13).  time  may 
be  computed  from  Eq  (3-16).  The  time  so  obtained  is 
normally  expressed  in  the  units  used  in  the  gravitational 
constant  ju  (See  Fig.  11). 

The  range  (radial  distance)  at  any  point  on  the  free- 
fall  trajectory  can  be  easily  found  after  one  point  of  the 
orbit's  eccentricity  circle  is  associated  with  a  range 
value.  Knowing  range  r,  at  point  (X,  ,Y,  )  the  square  of 
the  angular  momentum  per  unit  mass  is  obtained  from 


Eq  (3-10) 


h*  ~  y;  aa.  X, 


(3-17) 


The  second  point  fixes  Xt  and  by  the  conservation  of 
angular  momentum 
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Figure  11 .  Time  Lines  of  Constant  K 

(From  Ref  5*39) 


n  X2  -  ^ 


Thus  the  calculated  range  Is 
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It  follows  that  the  peripoint  radial  distance  must  be 


rX 
I  +■  e 


(3-20) 


The  semi-major  axis  for  all  conics  can  be  expressed  as 


a.  = 


rX 


I-  ez 


(3-2i) 


The  seml-latus  rectum,  p,  is  simply 


p  =  rX 


(3-22) 


Concept  and  Method  of  Graph  Usage 

In  order  to  increase  the  accuracy  of  data  extracted 
from  the  graphs  it  ie  evident  that  the  graph  must  be  quite 
large  and  have  a  denso  grid  of  lines  to  minimize  the  er¬ 
rors  associated  with  Interpolation.  Obviously,  combining 
cartesian  coordinates,  polar  coordinates,  and  time  con¬ 
stant  curves  on  one  graph  would  be  shear  nonsense. 

Instead,  it  is  suggested  that  there  be  three  separate 
graphs  plus  a  transparent  overlay  which  would  merely  serve 
to  transfer  the  plotted  points  from  one  graph  to  another. 

As  seen  in  Figure  12,  the  three  graphs  are  the  (X,Y) 
cartesian  graph,  the  (e,  Q  )  polar  graph  and  the  graph  of 
time  constant  curves  K.  If  the  distances  X  =  1  and  e  =  1 
are  represented  by  a  length  of  20  centimeters,  each  graph 


30 


GA/MC/71-1 


Transparent  Overlay 


lime  Constant  Curves,  K 


Figure  12,  Composite  Graphs 
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similar  to  those  shown  in  Figure  12  would  measure  50  cm 
horizontally  by  60  cm  vertically.  Line  Intervals  of  2  mm 
would  provide  X,Y  and  e  scale  readings  to  0.01  without 
Interpolation  and  to  about  0.001  with  interpolation.  In 
addition,  the  cartesian  and  polar  graphs  have  the  potential 
to  be  contracted  in  scale  by  factors  of  ten  in  order  to 
deal  with  orbits  having  an  eccentricity  less  than  0.1. 

Case  Studies  of  Applicability  to  Orbit  Determination 

For  simplicity,  let  it  be  assigned  that  a  set  of  three 
observations  exists  containing  range  r  (i.e.,  radial  dis¬ 
tance  from  the  center  of  the  attracting  body)  and  angular 
data  ^ where  (y  is  the  central  body  angle  or  great  circle 
arc  on  the  celestial  sphere  between  the  line-of-position  1 
and  the  succeeding  LOP's.  The  data  is  (r,  ,  %)  (tv  ,  ), 

and  (r,  ,  ) •  From  the  polar  coordinate  equation  for 

conic  sections,  one  can  derive  an  expression  for  X  and  Y 
in  terms  r*s  and  J^s.  (See  Appendix  E) .  With  X,Y  and  r 
known  the  problem  is  essentially  solved. 

A  second  case  can  be  considered  where  the  angular 
measurements  are  not  included.  Here  the  data  consists  of 
(r,  ,  t,  ) ,  (rt ,  tt )  and  (r3  ,  t3 )  where  t(  <  tz<  t^  .  Using 
equal  time  intervals  and  numerical  differentiation  (See 
Appendix  B) ,  it  is  possible  to  determine  r  and  r,  but  i.ot 
0  thus  posing  soma  difficulty.  The  difficulty  is  im¬ 
mediately  resolved  by  substituting  the  following  equation 
(See  Appendix  C) 
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jr©)  s  Vtc  +  v-ir  (3-23) 

into  Eq  (3-10)  producing  the  relationships 

x  =  i +  -5=?  (>24a) 


or 


J 


X 


i  +• 


rr 


(3-24b) 


where  is  the  local  circular  Telocity  for  radial  dis¬ 

tance  r.  Manipulation  of  Eq  (3-10) »  and  (3*11)  yields 


Y 


(3-25a) 


or 

(3*25b) 

The  solution  for  finding  the  geometric  orbital  parameters 
Immediately  follows. 

The  third  case  occurs  when  the  following  data  is 
obtained i  (  @,  »  t,  )  •  (Tj,  ,  »  ta )  and  (  t3  ) .  Here, 

X  can  be  found  from  Eq  (3-8)  and  from  the  first  time  deriv 
ative  of  the  angular  momentum  the  following  express lor  for 
r  is  given  as 
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(3-26) 


Having  found  r,  Y  can  be  solved  so  that  knowing  r,  X»  Y 
the  problem  is  solved. 

The  final  case  Is  the  situation  where  ranges  cannot 
be  obtained!  only  the  angular  measurements  are  available. 
The  data  consists  of  (  Qo  ,  t0  ) ,  (  0,  *  t,  ) ,  (  0,  *  t*)» 

(  Q3  ,  t3  ) ,  and  (  0+ ,  \  ) .  It  is  now  necessary  to  lean 
heavily  on  equations  derived  by  Hssrs.  R.H.  Gersten  and 
Z.E  Schwarz be In  (Hef  22i697).  Their  equations  arei 

e  cos  0(  =  0/S  (3-27) 

e  sU0,  =  R/5  (3-28) 


where 

Q  =  s t  n  ^a#3  —  (©3  4S m&0(* 

+  (@i  ]^[sm^0<2  -  sinA0f3j  (3-29) 

R  =  -(e/-^|cos60ra 

♦(s*-®/1)  tosae,3 

-(^-Ci  (3-30) 
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s  =  (®. sU (&e,.-ae.«) 

+  (ej*  su^e,t-(4)vSinAe„  o-%) 

Aft.  =  a  -ft  (3-32) 

‘j  * 


From  Eq  (3-27)  and  (3-28)  it  is  clear  that 


X  =  l  +  Q/s 


(3-33) 


and 

Y  =  R/5  (3-34) 

Although  X  and  Y  can  be  found,  a  range  must  be  associated 
with  the  point  (X,Y)  in  oilier  to  progress  with  the  solu¬ 
tion.  An  expression  for  range  can  be  obtained  from  Eq 
(3-8) »  it  is 


With  X,  Y,  and  the  corresponding  range,  r,  ,  known  the 
problem  is  readily  solvable. 


Summary 

The  main  advantage  of  the  velocity  hodograph  is  that 
by  mapping  all  conic  sections  into  circles  a  quick  graphi¬ 
cal  method  is  provided  for  the  determination  of  three 
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orbital  parameters  and  other  relationships.  The  orbital 
elements  are  the  eccentricity 

e-  read  directly  from  graph 
the  semi-major  axis 


a  = 


(3-21) 


and  time  to  or  from  peripoint 


At  =  K 


(3-16) 


Other  relationships  of  Interest  that  can  be  obtained  are 
the  true  anomaly 

0  -  read  directly  from  graph 
the  radial  distance  at  peripoint 

rj,  =  — X_  (3-20) 

r  i  +  e 

and  the  seml-latus  rectum 

p  -  T  X  (3-22) 

Other  advantages  are  that  the  hodograph  Is  Independ¬ 
ent  of  the  central  body  (l.e.,  It  can  be  used  when  orbit¬ 
ing  any  central  body,  hence,  serves  all  planets  and  the  sun) 
and  that  It  can  be  used  to  calculate  guidance  maneuvers 
(Ref  21«905)*  No  aspects  of  guidance  are  considered  in 
this  thesis  other  than  to  state  here  and  now  that  guidance 
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consists  of  orbital  maneuvers  based  on  conclusions  gained 
from  knowing  the  orbital  parameters.  Additional  Informa¬ 
tion  can  be  obtained  In  references  1  and  21. 

In  the  discussion  of  applicability  it  was  shown  that 
It  is  possible  to  obtain  from  range  and  angle,  range  only 
and  angle  only  data  the  information  required  to  fix  a 
point  (X,I)  on  the  hodograph  with  an  associated  value  of 
range. 
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IV.  Sample  Problems  and  Solutions 


Introduction 

Previously  in  Chapter  II,  the  position  fix  consisting 
of  a  line-of -position  (LOP)  and  a  radial  distance  was 
discussed.  The  LOP  is  a  line  from  the  attracting  body 
center  through  the  spacecraft  position  and  extended  to  the 
celestial  sphere.  The  intersection  of  the  LOP  with  the 
celestial  sphere  can  be  given  In  terms  of  Latitude  S  and 
East  Longitude  Since  the  radius  vector  contained  In 

the  orbital  plane  traces  an  arc  of  a  great  circle  on  the 
celestial  sphere,  two  orbital  orientation  parameters  — 
Inclination  and  right  ascension  — •  can  be  obtained  using 
spherical  trigonometry.  Furthermore,  knowing  the  true 
anomaly  which  can  be  found  from  the  velocity  hodograph,  in 
addition  to  knowing  the  inclination  and  right  ascension 
provides  the  necessary  data  to  determine  the  argument  of 
perlpoint  through  more  spherical  trigonometry  (See  Appen¬ 
dix  D).  The  three  orbital  parameters  —  inclination, 
right  ascension,  and  argument  of  perlpoint  —  will  not  be 
discussed  in  the  following  problems  since  the  calculations 
must  always  be  based  on  two  LOP's  and  follow  the  same  or 
comparable  mathematics  suggested  in  Appendix  D.  In  the 
sample  calculations  to  be  discussed  next,  a  cluster  of 
three  to  five  observations  is  typical  as  shown  in  Figure  13 . 
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The  central  body  in  all  cases  discussed  is  the  earth  where 

3  2 

the  gravitational  constant yx  is  62750.717  NM  /sec  . 

Sample  Problem  1_ 

Five  LOP's  have  been  established  at  15  minute  inter¬ 
vals  and  the  great  circle  arcs  separating  each  point  has 
been  calculated.  Since  angular  differences  are  involved 
let  if'.-O  degrees.  The  angular  data  is  #=  0,  /£=  6.464, 
12.787,  /£=  18.992.  and  %  =  25.106.  In  addition, 
one  range  sighting  was  obtained  on  the  third  observation) 
thus,  r3  -  15.475  NM 

From  the  numerical  differentiation  formulas  in 
Appendix  B,  the  time  derivatives  of  ^  expressed  in  sec¬ 
onds  are 


_  -2  5~./Q6  +  8  IZ.?8 

IZ  j  9 <30  ) 


i  69. 6  *  to  Jey/sec 
-  1.2 15  x  to*  Tra.d/sec 

and 

^  ..  1 1  ^  +  +  K}  -30  % 

IZ  t2 
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Plotting  X, 


~2f,/Q6  +  / 6  (zS.+Zt)  - 3o{tZ.787 ) 
/a  (ioo)(loo) 


-0.106  +  /6  (o.+se)  -3o(o.fcg7) 
12  (*Oo)(*OoJ 


=  -1.4-55'  *  !0  7  <fe§/sec 


2..S6  x  IQ  -ro.<J/seo 


(I.S+7S  x  I0+K-Z.S6  x  to') 
2(1,215  *  /O"4) 

-  /.6Z5  X  to"1  Nrt/sec. 

X.-rili 

(/■547g  x  /04j3(/.2/y  x  /Q~4^ 
6 . 275"  x  /04 

=  8.75  k  /o'1 


(3-26) 


(3-8) 


nr  -r*8 


(3-9) 


(1.625  <  t6')(i.5-+75  X  lO+fjl.ZIS  *  /o'4) 
6,275  #  /0  + 


=  7  56  *  /0’4 


Y  on  the  hodograph,  Fig.  14,  produces  an  eccen- 
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Figure  14,  Partial  Hodograph  with 

Sample  Problem  Solutions 


&  Problem  1 


□  Problem  3 


0  Problem  4 


O  Problem  2 
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tricity  of  about  .i5  and  a  true  anomaly  of  148  degrees. 
A  more  precise  graph  would  have  yielded 
e  =  0.1496 
0  =  147.5° 


a_ 


rX 

l-e* 


(3-21) 


U.S+7S  x  IO*)(8.7S  k  JO'') 
l  -(0./476J2 


-  J; s,  830  A/M 


The  exact  values  used  In  setting  up  the  problem  are 
a  =  13*760  NM 
e  =  0.15 
0  *  I5O.6960 


Sample  Problem  2 

Three  range  measurements  separated  by  15  minute  Inter' 
vals  are  obtained,  but  no  LOP  Is  established.  The  ranges 
expressed  In  nautical  miles  are  r,  =  9896,  rt  =  11098  and 
r3  =  12253.  Then  from  numerical  differentiation 


-  V 

at 


122.53  -  9896 
2  (90o) 

1. 307  Nh/sec 
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”  __  r-,  +  -  an 

*  ~  tz 


9896  +  IZZS3  -  2(11098) 
(loo)(loo) 


-  -S.8  *  /o"5  NM/*ecZ 


X=  I.+ 


(3 -24a) 


/  +  (/  <093  x  io*)z(-S.S  x  /o~5) 

'*  6.  £75*  X  /0+ 


=  0.886 


(3-25a) 


-  /.  309 


f(/./o9  8  >*  ro* )(0- 886) 
6-  275*  x  /04 


=  0.5*9 


Plotting  X,  Y  on  Figure  14  appears  to  yield  an  eccentricity 
of  0.526  and  102°  for  true  anomaly 


a  = 


(3-21) 
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I 


_  Q.foqe  *  io+)(o.s8b) 

I  -  0.526 a 

r  i3  6,oo  NM 

The  actual  values  used  in  this  problem  were 
a  =  13.760  NM 
e  =  0.5 
9  =  98.060° 


Sample  Problem  2. 

The  problem  to  be  now  considered  is  the  situation 
where  range  data  is  not  available  for  some  unexplained 
reason.  The  angular  data  obtained  over  20  minute  intervals 
is  =  0°,  0,=  12.883°,  25. 109°.  36.9440  and 

0^  =  48.657°.  The  angular  velocities  are  found  to  be 

0,  =  1.820  x  10-4  rad/sec 
# 

0t  =  1.745  x  10  rad/sec 
Qf  =  1.713  x  10"4  rad/sec 


Using  Eq  (3-29),  (3-30),  and  (3-31)  is  somewhat  laborious. 
In  order  to  retain  significant  digits,  a  desk  calculator 
was  used  to  accomplish  much  more  rapidly  what  could  have 
been  done  otherwise  by  hand  calculations.  The  results 
obtained  were 

Q  =  -2.5704  x  10“^ 
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E  =  1.5072  x  10'5 
3  =  1.4939  x  10”4 
and  from  Eq  (3-33)  and  (3-34) 

X  =  .8279 


Y  =  .1009 


The  range  corresponding  to  position  0 t  is 


0-35) 


and  yields  a  value  for  r,  of  11575  NM. 

Vlhen  X,Y  is  plotted  on  Figure  14  it  appears  to  indi¬ 
cate  an  eccentricity  of  0.2  and  a  true  anomaly  of  150°. 
The  actual  values  in  setting  up  this  problem  are 


e  —  0.2 
0(«  150.658° 
r(  =  11,62?  NM 
a  =  10,000  NM 


Sample  Problem  4 

This  last  example  represents  the  case  where  three 
complete  position  fixes  are  known.  The  data  has  been 
reduoed  to 
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11,489  NM 

tl 

0 

0 

12,604  NM 

^  =  7.541° 

13,619  NM 

^  =  14.482' 

The  equations  to  be  used  are 


- 

r*.  -  X 

— 

Ti  -  17 

_  77  -  cos  ^ 

"T,  -  -r3  cos  % 

*  tin 

- 

13  «*«-  *6 

TT  ~  ij  cos 

T  -  ’Jco./j 

(E-8) 


and 


e  cos  Qt 


*»•  -  *r 

sin  ^ 

.’r 

[T7-1JCOS^_ 

sln.Qt 


(K-9) 


See  Appendix  E  for  derivation  of  the  above  formulas. 
Breaking  down  Eq  (E-8)  into  its  basic  terms  yields 


ZL-Z-2L  —  /2.fr04  -  //489 

17  -^Cos/%  //489  -  (i2-60+)(o.Wf3S) 
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ill  6 
-  I  o©  S 


I.IO 9 


n  -  y,  _  *36(9-  tf±si 

X  -X  cds  ft*  //+ 8?  -(t36ll)(o,96823) 

I  4  J 


2130 
-  J707 


-  ~l.Z+f 


rasin  /fr,  _  0 z<,o+)(o./3iz. *) 

X-Tz  cos  ^  "  —loo  S’ 


=  -A  65*0 


Tg  S6W.  ^  _  (/36/?Xo*2gQOg) 

V-TJCOS^  "  ”  1707 


-  -1.990 
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e&inft  £ 

1  -f.GSO  +  f.990 


0.4-00 


ecos©,  -  -4/04  -(-/.fes'o )(o.+oo) 


-  -o.  4-4*9 


Therefore 


X  =  0.551 


o 


and 


Y  =  0.400 


Plotting  X.Y  on  Figure  14  produces  an  eccentricity  of 
about  0.6  and  a  true  anomaly  of  about  138  degrees. 

The  semi-major  axis  is 


Gt  = 


(3-21) 


Ql+8l)(0S5-l) 

I  -  (o.t»)z 


/  0,000  NM 
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The  actual  values  in  this  problem  are 

e  =  0.6 
0  =  137.580° 
a  *  10,000  NM 


Discussion 

The  sample  problems  were  Intended  to  show  the  reader 
that  situations  where  the  data  is  restricted  In  some  form 
or  other  can  be  solved  Just  as  the  theory  Indicates.  Ad¬ 
mittedly,  the  examples  selected  to  demonstrate  the  various 
situations  were  chosen  with  care. 

Sample  problems  1,  2,  and  3  rely  on  derivatives  ob¬ 
tained  through  numerical  differentiation  and  consequently 
contain  truncation  errors.  Truncation  errors  are  mini¬ 
mized  as  the  time  Interval  gets  smaller.  On  the  other 
hand,  measurement  errors  are  magnified  by  the  numerical 
differentiation  algorithms  as  the  time  interval  gets 
smaller  (Ref  17*9'*'98).  All  the  sample  cases  selected 
consisted  of  range  data  correct  to  the  nearest  nautical 
mile  and  angular  data  correct  to  the  nearest  thousandth  of 
a  degree.  Experimentation  with  time  intervals  of  5.  10, 
15,  20,  and  30  minutes  and  also  experimentation  with  3  and 
5  point  differentiation  did  not  vary  the  results  of  eccen¬ 
tricity  and  semi-major  axis  by  more  than  2  per  cent  when 
a  fixed  position  on  the  orbit  was  selected.  By  contrast, 
however,  the  point  of  the  orbit  selected  for  experimenta- 
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tion  had  an  enormous  effect.  Had  Sample  Problem  2  been 
set  up  with  the  same  time  interval  but  at  the  true  anomaly 
of  150°  instead  of  98°*  the  solution  would  have  produced 
an  eccentricity  of  about  0.4  and  a  semi-major  axis  of 
about  15.500  NM  whereas  the  true  values  were  0.5  and 
13.760  NM,  respectively.  Obviously  this  represents  a 
deviation  which  cannot  be  neglected,  not  even  by  the  most 
careless  navigator.  That  approximate  derivatives  obtained 
from  finite  differences  should  be  viewed  with  skepticism 
(Ref  17«97)  seems  to  be  good  advice  until  some  future 
study  reveals  the  criteria  by  which  the  e*.d  results  ob¬ 
tained  through  calculations  involving  numerical  differen¬ 
tiation  can  be  accepted  or  rejected. 

Sample  Problem  4  consisting  of  three  complete  position 
fixes  does  not  rely  on  derivatives  and  not  surprisingly 
should  yield  an  acceptable  result  always. 

In  Sample  Problems  1  and  3.  the  shortcomings  of  try¬ 
ing  to  retain  significant  digits  through  slide  rule  cal¬ 
culations  were  quite  evident  to  this  author.  To  have 
measuring  instruments  capable  of  measuring  angles  to  the 
thousandth  of  a  degree  and  using  five  place  trigonometric 
tables  but  not  being  able  to  retain  significant  digits  in 
the  calculations  seems  somewhat  unreasonable.  A  hand-held 
mechanical  calculator  such  as  the  Curta  Calculator  which 
is  smaller  than  a  four  inch  cube  and  weights  lass  than  one 
pound  yet  has  an  eight  significant  figure  accuracy  (Ref  61 
167)  appears  to  be  a  necessary  item  for  any  mission 
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V.  Gro ttrd  Based  Orbit  Determination 

Introduction 

Since  electrical  navigational  devices  are  barred  from 

consideration  in  this  thesis,  the  tracking  of  satellites 

VV. 

must  depend  on  optical  means.  In  order  to  be  visible,  the 
spacecraft  must  carry  Its  own  light  source  or  be  Illumi¬ 
nated  by  an  external  source  such  as  the  sun. 

The  optical  sightings  normally  consist  of  angular 
measurements  In  azimuth  and  elevation  taken  from  a  ’"sen¬ 
sor"  site  local  horizon. 

Certain  assumptions  to  be  made  are  that 

1.  The  coordinates  of  the  site  location  are 
known. 

2.  The  direction  of  planetocentric  true  "North" 
is  known. 

3.  The  true  local  vertical  can  be  established, 

4.  The  central  body  surface  radius  is  known. 

5.  The  central  body  rotation  rate  is  known. 

Simple  Devices 

The  devices  mentioned  are  crude  devices  in  most 
instances  and  to  a  great  extent  resemble  the  schemes  used 
by  amateur  moonwatchers . 
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The  Mast 

The  mast  Is  the  simplest  device.  It  consists  of 
erecting  a  mast  of  known  height  with  its  longitudinal  axis 
parallel  to  the  local  vertical.  The  ground  about  the  site 
should  be  level  to  approximate  a  plane  perpendicular  to 
the  local  vertical. 

The  individual  making  the  sightings  must  move  about 
the  mast  in  order  to  maintain  the  visible  satellite’s 
position  just  on  top  of  the  mast.  As  the  Individual 
tracks  the  satellite  by  moving  his  position  about  the 
mast,  marking  flags  are  dropped  at  recorded  times.  After 
satellite  passage,  the  dropped  flags  can  be  used  to  di¬ 
rectly  establish  azimuth  and  time.  Knowing  the  observer’s 
eyesight  height  and  a  little  plane  trigonometry  can  estab¬ 
lish  the  angle  of  elevation  (Ref  7:62).  See  Figure  15. 

Slew  Telescope 

The  slew  telescope  is  essentially  the  most  common 
instrument  used  to  simultaneously  provide  a  reading  of 
azimuth  and  elevation.  The  principle  is  quite  simple, 
featuring  two  degrees  of  freedom.  In  essence,  the  instru¬ 
ment  consists  of  a  sighting  devioe  mounted  on  a  protractor 
which  is  free  to  swing  cn  a  pin  connected  to  a  vertical 
rotary  shaft.  A  sketch  Is  illustrated  in  Figure  16. 

Techniques 


In  order  to  avoid  using  a  laborious,  classical  orbit 
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Figure  15.  The  Mast  Figure  16.  Slew 

Telescope 

determination  scheme  such  as  Laplace's  Method,  the  slant 
range  to  the  target  must  be  found  to  simplify  calculations. 
The  methods  used  to  find  the  range  resort  to  triangulation, 
but  all  triangulations  require  one  known  distance  called 
the  baseline. 

There  are  two  basic  techniques  available.  They  are* 

1.  The  two  station-one  pass  method 

2.  The  one  station- two  pass  method. 

The  two  station-one  pass  method  can  establish  the 
range  to  a  satellite  from  only  one  satellite  passage.  The 
baseline  is  the  known  separation  distance  between  the  two 
stations.  Due  to  restricted  mobility,  a  mission  landing 
on  some  planet  would  not  generally  be  expected  to  es tab- 
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Figure  I?.  Triangulation  Error 
(From  Hef  11:121) 

llsh  a  long  baseline.  If  the  baseline  Is  small  compared 
to  the  height  of  the  target,  the  closing  angle  between  the 
two  rays  Is  extremely  small,  and  slight  tracking  or  mis¬ 
alignment  errors  in  the  direction  to  the  target  will  pro¬ 
duce  tremendous  errors  in  range  and  spatial  position  to 
the  target  (Ref  11:119-120).  See  Figure  17. 

The  one  station-two  pass  method  makes  use  of  the 
fact  that  an  unperturbed  orbital  plane  remains  stationary 
but  due  to  the  planet's  rotation  the  sensor  site  has  moved 
Its  location  with  respect  to  that  plane  during  the  elapsed 
time  between  satellite  passage.  The  baseline  is  the  chord 
which  separates  the  site  locations  with  respect  to  the 
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orbital  plane  during  each  passage.  Obviously  the  distance 
of  the  baseline  is  a  function  of  angular  rate  of  planet 
rotation  and  elapsed  time.  The  time  interval  must  corre¬ 
spond  to  the  time  between  sighting  the  satellite  on  the 
first  pass  and  sighting  the  satellite  at  the  same  position 
in  its  orbit  on  a  succeeding  pass.  The  baseline  distance 
B  can  be  expressed  as 

B=  2Rkcosi  sin^  (5-1) 


where 

R  =  body  surface  radius 
S  =  site  latitude 

=  angular  velocity  of  planet  rotation 
=  elapsed  time 

From  the  spherical  isos¬ 
celes  triangle  of  Figure  18 
it  can  Ke  seen  that  the  angles 
A,  and  (36O0  -  Aa)  must  be 
equal  and  their  magnitude  is 
a  function  of  Atf.  With  a 
precomputed  table  containing 
A,  and  At  versus  values  of  &t, 
a  search  of  the  tracking  data 
can  be  made  to  find  the  azi¬ 
muths  A,  and  A.  which  will  Figure  18.  Azimuths 

Versus 

Elapsed  Time 
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correspond,  to  the  elapsed  time.  The  correct  correspond¬ 
ence  establishes  a  point  on  the  orbit  that  is  the  common 
reference  point  between  the  two  passages  of  the  satellite. 

The  greatest  difficulty  arising  from  triangulation 
with  the  baseline  is  that  the  angle  between  the  baseline 
and  the  satellite  line-of-slght  must  be  calculated  through 
vector  dot  products.  See  Appendix  F. 

Satellite  Position  Fix 

After  calculating  three  slant  ranges  by  triangala- 
tion,  the  satellite  position  can  be  determined  quickly. 

The  radial  distance  from  body  center  to  spacecraft  is 

T  -  ^  sin(El)  (5-2) 


where 

=  body  surface  radius 
rs  =  slant  range 
El  =  elevation  angle 

The  great  circle  arc  or  angular  ground  distance  4>  from 
the  site  to  the  sub-satellite  point  is  found  from 


*  *< 


-I 


an 


cos 


L  +  S  sin  (El) 


(5-3) 
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Figure  19  Figure  20 

Range-Elevation  Chart  Geometry  for 

Sub-Satellite  Point 

An  all-purpose  non-dimensional  chart  such  as  illustrated 
in  Figure  19  can  be  used  to  obtain  radial  distance  r  and 
angular  ground  distance  (p. 

Application  of  the  sine  and  cosine  laws  of  spherical 
trigonometry  allows  one  to  calculate  the  position  of  the 
sub-satellite  point  in  terms  of  latitude  and  longitude 
(See  figure  20).  Knowing  the  sub-satellite  point  is  anal 
ogous  to  knowing  the  LOP  in  space  navigation.  From  this 
point  on  the  solution  is  arrived  at  in  the  same  manner  as 
Sample  Problem  4  where  three  complete  position  fixes  are 
known. 
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Stannary 

The  foregoing  discussion  merely  served  to  show  that 
satellite  tiacking  does  not  differ  greatly  from  space 
navigation.  The  tricky  part  concerns  manipulation  of  the 
angular  data  to  establish  the  slant  range  and  then  to 
determine  the  sub-satellite  points.  Once  three  satellite 
position  fixes  are  known*  the  problem  is  solved  in  the 
same  way  as  for  space  navigation  where  three  complete 
position  fixes  are  known. 
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VI.  Conclusions  and  Recommendations 

Introduction 

This  chapter  comments  on  some  of  the  Impressions 
drawn  from  the  i*esearch  associated  with  this  thesis. 

Under  the  subheading  "Recommendations,"  some  discussion 
of  follow-on  theses  are  briefly  touched  upon. 

Conclusions 

The  following  conclusions  are  offered  as  the  result 
of  experience  gained  in  preparing  this  thesis* 

1.  Navigational  devices  which  are  simple  in 
concept  and  design  appear  to  present  no 
problem  from  the  theoretical  point  of  view. 
The  real  problem  arises  when  the  instrument 
is  Integrated  into  a  confined  system  which 
places  restrictions  on  size,  weight,  and 
employment.  Then,  the  simple  device  is 
transformed  into  a  complex  instrument  re¬ 
quiring  advanced  and  sophisticated  engineer¬ 
ing  deslgnj  witness  the  transformation  of 
the  simple  marine  sextant  to  the  complex 
space  sextant. 

2.  The  effects  of  oblateness  on  the  range 
accuracy  of  stadimetric  measurements  involv- 
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ing  devices  viewing  the  entire  target  can  be 
resolved  by  additional  calculations.  Thus, 
the  accuracy  in  range  is  really  dependent  on 
how  well  the  physical  dimensions  of  a  celes¬ 
tial  body  are  known. 

3.  The  composite  graphs  of  the  non-dimensional 
velocity  hodograph  appear  to  provide  the 
navigator  with  a  tool  which  can  assist  him 
in  determining  three  orbital  parameters 
(eccentricity,  semi-major  axis,  and  time  to 
or  from  peripoint  passage)  plus  the  true 
anomaly.  In  addition  to  the  fact  that  the 
velocity  hodograph  is  a  simple  picture  of 
any  conical  orbit,  the  hodograph  has  the 
potential  of  being  used  for  guidance  or 
orbital  maneuvers  which  are  usually  the 
normal  follow-up  actions  after  determina¬ 
tion. 

4.  Orbit  determination  can  be  completely  or 
partially  accomplished  depending  on  what 
data  is  available  by  using  first  and  second 
time  derivatives  of  range  and  true  anomaly 
obtained  through  numerical  differentiation. 

5.  The  accuracy  of  the  geometric  orbital  para¬ 
meters  derived  through  numerical  differen¬ 
tiation  is  more  sensitive  to  the  position 
of  the  spacecraft  on  its  orbit  than  to  the 
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time  interval  chosen. 

6.  Manual  navigation  aids  such  as  a  mechanical 
calculator  is  a  necessity  in  order  to  retain 
the  significant  digits  within  the  computa¬ 
tions. 

7.  Ground  level  orbit  determination  is  basical¬ 
ly  no  different  from  on-board  orbit  determi¬ 
nation.  Without  resorting  to  Laplace's 
Method  or  some  similar  scheme,  the  key  to 
the  problem’s  solution  is  determining  range 
by  triangulation  so  that  three  spatial  fixes 
are  obtained. 


Recommendations 

Some  suggestions  are  offered  for  follow-on  study 
topics.  They  are  as  follows i 

1 .  Resolve  the  geometry  of  the  stadimetric 
system  viewing  a  limited  portion  of  a  non- 
spherical  target. 

2.  Develop  large  size  composite  graphs  of  the 
velocity  hodograph  including  a  finely  in¬ 
cremented  constant  time  line  graph  and  to 
investigate  the  handiness  of  the  composite 
graphs  thoroughly. 

3.  Investigate  what  criteria  should  be  used  to 
accept  or  reject  the  geometric  orbital  para¬ 
meters  obtained  from  calculations  relying 
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Appendix  A 


Supplementary  Relationships  in  the  Geometry 
of  Stadlmetric  Measurement 


The  equations  below  (Ref  13*101)  refer  to  Figure  5. 
The  angle  Q-&Q  is  measured;  &Q  is  the  correction  angle 
due  to  the  finite  size  of  the  center  field  of  view. 


Let  R  =  r^,  then 


~y  -  SlU 


iK 

'"( i+Vfl) 
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*  s  7 
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=  COS 


( $>  2  A 
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Appendix  B 

Numerical  Differentiation  Formulas 

The  first  and  second  time  derivatives  are  given  below 
for  finite  differences  obtained  from  three  to  five  read** 
ings. 

3  readings i 


*3  -  *, 
2 1 


(B-l ) 


*  • 


X,  +  X3  “  2 
t2 


(B-2) 


4  readings t 


3  (  ^  *3  ~  +  *♦) 

6t 


(B-3) 


X,  +  X3  -  2  X. 

tz 


(B-4) 
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5  readlngst 


(*,-*«)  -  8  (xt  -X4) 

IZ  t 


(B-5) 


v  _  -(\  +  xJ-F/«fx4+Xt)-30Xj 

a  -  j^jz  (B-6) 
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Appendix  C 


Radial  Acceleration 


When  a  body  is  in  an  exact  circular  motion  it  may  be 
said  that  the  centripetal  acceleration  is  balanced  by  the 
gravitational  acceleration  in  which  case  the  radius  is 
constant  implying  that  r  =  r  =  0.  Thus 


(C-la) 


or 


(C-lb) 


where  g  is  the  local  gravitational  acceleration  at  dis¬ 
tance  r 

ls  the  local  circular  angular  velocity 

VLC  Is  the  local  circular  velocity 
•  -  • 

When  e*etc  ,  the  accelerrtions  are  not  in  balance » 
and,  there  exists  a  radial  acceleration  r  and  consequently 
a  radial  velocity  r.  The  velocity  is  expressed  as 
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u 

and  the  relationship  for 

radial  acceleration  as 

r0 

l 

•  • 

(C-3a) 

or 

f  • 

z 

•  • 

[re 

3-^  +  rir 

(C-3b) 

Substituting  Eq 

(C-lb)  into  (C-3b)  yields 

2 
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H)  +  rlr 

(C-4a) 
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Appendix  D 


Orbital  Orientation  Parameters 


It  is  assumed  that  two 
lines  of  position  have  been 
identified  in  terms  of *  and  i. 
The  most  general  equations  for 
the  orientation  parameters 
(Hef  10 1265)  are  given  below. 

Right  ascension, -A.  i 


Figure  21 
Three  Dimensional 
Geometry  of  an  Orbit 


f 

-fl.  =  tou\  '  { 


tan  st not. 

Uni, 


cos*,  -  cos*2 
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Inclination,  ii 


(D-2) 


The  true  anomaly  angle  between  the  reference  plane 


and  point  1  is  given  by 
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.  -i 

sin. 


sin.  &, 
sin.  i 


from  which  the  argument  of  peripoint  is 


60  - 


where  is  the  true  anomaly  at  point  1  as  obtained, 
the  velocity  hodograph. 


(D-3) 


(D-4) 

from 
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Appendix  E 

Geometric  Parameters  from  Range 
and  Angle  Data 

This  method  can  be  used  to  find  X  and  Y  for  the  ve¬ 
locity  hodograph  when  three  complete  position  fixes  are 
known*  Since  the  LOP*s  of  the  position  fixes  are  expressed 
in  terms  of  latitude  &  and  longitude <*,  the  great  circle 
arcs  connecting  LOP's  with  the  LOP  of  position  fix  1  form 
spherical  triangles  where  the  two  sides  (90°  -  and  90° 

-  /.  )  and  the  included  polar  angle  (  -  o<t  )  are  all 

known.  As  defined  here,  the  great  circle  arc  of  LOP  1  is 
a  special  case  and  may  simply  be  expressed  as 

p  -  O  (E-l ) 

The  succeeding  great  circle  arcs  may  be  found  from  the  law 
of  cosines,  thus 


^  =  cos"7  [  cos  (10-St)  cos(w  '«fi) 

+  Sin  (e-2) 

which  can  be  simplified  to 

s  cos  1  [sin.J,  fiU/j  +  cos^cos^  cos(«<.  - «Ce)]  (e-3) 
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It  is  evident  that  the  time  separation  between  fixes  is 
immaterial  in  the  above. 

From  the  general  polar  equation  of  a  conic  trajectory 
Y"(/+ecos0 )  =  constant  (E-4) 

therefore  let 

r  (l  +  e  cos0j  cos(0  +  £)J  (E“5) 

and 

^  (i  +  e  cos  6)  =  ir3(i  +  e  cos(0  +  f>))  <E“6) 

Solving  simultaneously  by  making  use  of  the  trigonometric 
identity 

COS  (9-*-  a  COS0COS^  -sinQsinft'  (E-7) 


yields 
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and 


e  COS  Qt 


r  -  v 

I"  TZsinPL  1 

|jr-T£cos^_ 

”  W-ricos  ^  J 

e  sin  9, 


(E-9) 


Although  Eq  (E-8)  and  (E-9)  look  long,  the  calculations  In 
Eq  (E-8)  are  simple  and  straightforward*  and,  the  terms  of 
Eq  (E-9)  are  repetitions  of  terms  In  Eq  (E-8).  Note  that 

=  /  +  e  cos9t  (3-10) 

and 

Y  r  e  sin-0,  (3-11) 

are  the  arguments  for  the  velocity  hodograph. 
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Appendix  P 

Baseline  to  Llne-of-Slght  Angle  ( 3~) 

To  find  the  slant  range  along  the  llne-of-slght  one 
must  know  the  angle  formed  between  the  baseline  and  llne- 
of-slght  before  trlangulatlon  car  be  made  with  the  base¬ 
line  distance.  A  coordinate  system  Is  established  such 
as  Illustrated  In  Figure  22. 

The  llne-of-slght  unit  vector  can  be  expressed  as 

Ls  r  Cos(A*)cos (zi)  X -sm(As)cos(El)ij  +  sin(£i )j,  (f-1  ) 


Figure  22 

Llne-of-Slght  Unit  Vector 


( 
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